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Abstract 

In this paper we consider generalized eigenvalue problems for a family 
of operators with a quadratic dependence on a complex parameter. Our 
model is L(A) = -A + {P{x) - \f in L^(R'') where P is a positive elliptic 
polynomial in R'' of degree m > 2. It is known that for d even, or d = 1, 
or d = 3 and m > 6, there exist A G C and u £ L'^(R'^), it 7^ 0, such 
that L{X)u = 0. In this paper, we give a method to prove existence of 
non trivial solutions for the equation L{X)u = 0, valid in every dimension. 
This is a partial answer to a conjecture in ^12^. 

1 Introduction 

Let us introduce the following family of differential operators, 

L{X)^-A, + {P{x)^Xf (1.1) 

where A^; is the Laplace operator in R^, A is a complex parameter, P is a 
polynomial of degree m > 2 such that the leading homogeneous part P,„ of P 
satisfies Pm{x) > for every x E R''\{0} (in other words we say that P is a 
positive-elliptic polynomial) . 

Such family of operators play an important role when studying analytic smooth- 
ness of solutions of differential operators with multiple characteristics (see [T^ 
and references there). They also appear in the theory of damped oscillations 
in mechanics [21 [H]. The question we want to adress here is: "does there exist 
A G C and u in the Schwartz space iS(K''), u ^ 0, such that L{X)u = 0?" In 
[4] the authors have proven existence of non trivial solutions for 1 < d < 3, 
assuming that m is large enough for d = 3. After, Helffer-Robert-Wang proved 
in [12] the following result. 
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Theorem 1.1 Assume that d is even and that P is a positive- elliptic polyno- 
mial of degree m > 2. 

Then there exist A G C and u G iS(M'^'), w 7^ 0, such that L{X)u = 0. 

The proof given in [12j shows that there exist an infinite number of such eigen- 
values [in] located in the half-plane {A G C, 5RA > 0}. But it is not known if 
the generalized eigenfunctions span all the Hilbert space i^(R'^), excepted for 

d = 1 [ig. 

For d odd , d > 3, m > 2, the problem of existence of non zero solutions is 
still open and it was conjecture in T5] that such solutions exist whatever the 
dimension d. 

In this paper we prove that this is true for every elliptic polynomial if d = 3 and 
for large classes of elliptic polynomials for d = 5,7. We also discuss a numerical 
approach to prove that some coefficient in a semi-classical trace formula is not 
zero. For d > 9 we conjecture that this coefficient is not zero hence there exist 
an infinite number of nonlinear eigenvalues. 

This work was supported by the program ANR 08-BLAN-0228, NONAa, Re- 
search French Ministry. 



2 nonlinear eigenvalue problems 

In this section we recall some known properties concerning nonlinear eigenvalue 
problems. For more details we refer to pUl fT51 [TO] , 

Let us consider the quadratic family of operators L{X) = Lq + XLi + A^ where 

Lq, Li are operators in an Hilbert space Ti.. Lq is assumed to be self-adjoint, 

1/2 

positive, with a domain D{Lo) and Li is vXo-bounded. Moreover Lq is in 

a Schatten class C^iji.) for some real p > 0. 

The following results are well known. 

Theorem 2.1 L{\) is a family of closed operators in Ti.. 
A I— > i^^(A) is meromorphic in the complex plane. 

The poles Xj of L~^(X), with multiplicity ruj, coincide with the eigenvalues 
with the same multiplicities, of the matrix operator Al in the Hilbert space 
H X D{Ly^), with domain D{Al) — L){Lq) x D{lI^^) where 

A. - ( X 1 ) ■ <-) 



Let us denote Sp[L] the eigenvalues of Al (which coincide with the poles of 
L~\z)). 

Remark 2.2 It may happens that Sp[L] is empty. The following one dimen- 
sional example is interesting and was discussed in \16l O Idjj. 

L^.g{X)=-^ + {x'^ -Xf +gx"'-\ (2.3) 
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For every m>2,m even, Lm,o has infinity many eigenvalues but L^.m has no 
eigenvalue. The last statement is a conseguence of the factorization 

ax ax 

So, we can compute all solutions for the eguation Lm.m{^)u — and see that a 
non-null solution u is never hounded on M. 

But if m is odd, Lm,m{^)u = 0, has infinity many eigenvalues on the imaginary 
axis 

On the other side there exist sufficient general conditions to have Sp[L] ^ 
ll(K\15^ . Unfortunately these conditions are not fulfilled for our example L{\) = 
-A^ + {F{x) - \f when d>2. 

The following formula appears for the first time in [2] and will be very useful 
for our purpose. 

Theorem 2.3 For k large enough (k ^ N, k > p) and for z G C\Sp[L],, we 
have 

Tr{AL - z)-'^-' = ^Tii^{L{z)-'L'{z)], (2.4) 
where each above operators are trace class. 
Using Lidskii Theorem _10J and (|2.4p . we get 

Y: m.(A-z)-'=-i^^Tr[^(L(z)-iL'(z)))]. (2.5) 

AeSp[L] 

where m(A) is the multiplicity of the eigenvalue A. 

As it was nicely remarked in the paper |4j , a sufficient condition for Sp[L] ^ 
is that the r.h.s in (|2.5p is not zero. To check this property a natural method is 
to introduce parameters and use semiclassical analysis. 

In |T5] the authors also use Lidskii theorem and semi-classical analysis on the 
matrix system Al ■ Here we consider more directly the scalar family of operators 
L{z) where computations are easier even if the dependence in z is nonlinear. 



3 Semiclassical parametrix 

For simplicity we assume here that F is homogeneous of degree m > 2 and 
F{x) > for X e E^, X 7^ 0. By the scaling transformation with 

fi = ^-{m+i)/m g^^^ 2: — — we can see that L{\) is unitary equivalent to the 

T 

semiclassical Hamiltonian r^L(z) where 

L{z) = -h^A^ + {F{x) - zf . (3.6) 

L{z) is the ft-Weyl operator with the symbol L{z,x,£^) = + (F{x) — z)^. For 
semiclassical analysis tools and ?i-Weyl quantization we refer to \18^ . Here we use 



4 



the notation H for the S-Weyl quantization of the symbol H or for convenience, 
H = Opl{H). 

Using semiclassical operator calculus, we can construct a good parametrix 
for L{z)^^ for z e A where A is the sector 

A = {z e C, \z\ > To, tt/2 + S < arg(z) < 37r/2 - 5}; tq > 0, (5 > 0. 



Theorem 3.1 There exists a semiclassical symbol K^^\z), z £ A, < h < 1, 
such that 

L{z)-' ^ Opl{Kn{z)). (3.7) 

Moreover the asymptotic expansion has the following meaning: for every N > 1 
we have 

L{z).Op^l J2 h''K,,{z)\ =l + h'^+'Of^{4''^{z)) 

\0<j<N J 

where the symbol r!^2n{^) satisfies the following estimates : 
for every a, /3 G N'' we have 

< c(iv...ffl ^"-;^;;';;.!t^:;;^'"V (^.':)-"-'-'-'" 

(3.8) 

for every a,/? G N'', where C{N,a, f3) is uniform in z £ A and where fj,{x,^) = 
(l + |a:|2'" + |$|2)i/2m^ 

sketch of proof. The method to get such result is standard and was used 
many times to construct parametrix of elliptic pseudo-differential operators |20| . 
Usually the z-dependence is linear but here it is quadratic. Moreover here 
we need accurate estimates for the remainder term in the product of pseudo- 
differential operators depending on parameters. The necessary estimates for 
i?2^(^;^'0 ^i'^ established using the technics coming from the papers [71 [3]. 
An other difficulty here is that we shall need to compute the symbols K2j for j 
large enough. This computations are not easy, so we have to be explicite as far 
as possible. 

Using the product formula for ft-pseudodifferential operators, we get at the 
initial step: 

^ H^) ^ I^P + (4) - ^'-'^ 

and the induction formula 

K2,^-Ko[ J2 E na,P)d^d(^L{z)d^^d^K2i\ (3.10) 
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where T{a,l3) = 2^u-')a\i3\ ' compute K2 and K4. 

^ L3(z) L4(z)' 

L2{z) = {P{x)- z)AP{x) + \S/P{x)\\ 

L^iz) = -2[{P{x)-z)D'P{x)^-^+{VP{x)-0' + {P{x)~zf\VP{x)\% 

where D^P{x) is the Hessian matrix of P in variable x. 
Now using (|3.10p we have 

= -A'o{E|/3|=4r(o,/3)9^£(z)afi^o + E|a|=2r(a,o)a^"L(z)a,"if2 

+ E|/3|=2r(0,/3)9.^i(z)9fif2)}. (3.11) 
By induction on j, we easily get that 

j+l<k<3] ^ ' 

Q^J {x^ P — z,^) is a polynomial in ((P — z), ^), with a total degree < fc — 2, with 
coefRcients depending on derivatives of P{x). 

The following lemma will be useful later. Let us denote val[(5^"'], the valuation of 
Q^"* as a polynomial in P — z, ^. Let us recall the definition of valuation. Denote 
by I the ideal with generators <^i, • • • ,£,d,P — z, in the ring C°°{M.^ x M.^)- If 
Q e C°°{R^ xR^), val[Q] is the biggest integer p such that Q elP. 

Lemma 3.2 We have 

val[gf ] > 2(fc - 1 - 2j), for 2j + 2 < fc < 3j, and j > 1. 

Proof. This is easily proved by induction on j, using (|3.12p and the following 
formula. Let Q and L be smooth functions in R", a multiindex a G N", then 
we have 

(v^J"-^ ^ 

where in the sum we have the conditions, 7j G N", /ij G N, 7 < a, • • = /i, 
M1I71I H ^ ^^|>\l^\ = □ 

Remark 3.3 The parametrix computed above is enough to get qualitative infor- 
mations. Quantitative informations are much more difficult to get except for the 
first orders (j = When j is larger it is not so easy to compute explicitely 

the terms Q^f^ {x, P — z,£^). 
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Remark 3.4 It is not difficult to extend the above results when the elliptic 
polynomial P{x) has lower terms: P ~ Pm + Pm~i + • • • Pi + Po where Pj is 
homogeous with degree j and Pm{x) > for x € M'*\{0}. Then we have 

P(ri/™y) = rP(^)(y) 

with e = T-i/'" = and P^'Hy) = PUy) + ePm-i{v) + ■■■+ e^PoCy). 

So P(^^ is a uniform elliptic family of polynomials and we can easily see that 
the constructions in US. 8\) are uniform in the small parameter e. 



4 A trace formula 

Recall that Sp[L] denote the generalized eigenvalues of the quadratic family 
L{z), m\ is the multiplicity of the eigenvalue A. Let / an holomorphic function 
in A such that 

|/(z)| <C(l + |z|)-^ VzeA. (4.14) 

For our applications we shall choose f{z) = (z + A)^'', for a suitable parameter 
A e C. Let be F a complex contour in A defined as follows. 

F = {re±'«°, r > ro} U {roc^^ 0o < ^ < 2^ - 0o}, 

where ro > and ^ < 9q < tt. 

Proposition 4.1 Assume that /i > ^^^^^^7^- Then fiAh) is a trace class oper- 
ator and we have 

Tr(f(A)= mAf(A) =Tr[/L(z)-iL'(z)f(z)dz], (4.15) 

AeSp[L] ''^ 

where §^F{z)dz = J^F{z)dz (contour integral in the complex plane). 

Proof. This a direct consequence of the Cauchy integral formula and Theorem 
[231 □ 



Theorem 4.2 For f as above, for every d > 1 we have in the semiclassic regime 
h\0, modulo 0{h+°°), 

^ m,/(A)x^cg)(/)n2,-d_ 

AeSp[L] j>0 

If d is odd, 

d,'\f) = (4.17) 

and for d even, 

&^\f) = 2{~lY'\2n)-^ [ [ f{P{x) + |r;|)dxdr;. (4.18) 
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For the other terms (j > 1) we have the following qualitative information 

^2?(/)= E / ^2,,fe(x)/W(P(x))dx, (4.19) 

where A2j^k{x) are polynomials in d2P{x), I7I < 2j and nj depends on j . 



Moreover if d is odd, then c'^\f) — for d > 4j + 1 



Proof. The asymptotic expansion (|4.16p is a direct consequence of (|3.8[) and 
usual properties of trace operation for Weyl quantization. 
Let us compute Cq'^^(/). Wc have the integral formula: 

where dx — (27r)"'^o?x. By the residue theorem we get 

C^'^(/)= / / [f{P{x)+^\S,\) + f{P{x)^^md^dx. 
For a > we have 

(/(P(x) + a\imdx = a-^ [ [ {f{P{x) + l^Mdx. 



So by analytic extension and evalution at a = i we get formula (|4.17p and (|4.18p . 
In particular we see that for d even, there exists / satisfying (|4.14p such that 

For j > 1, using (|3.12p . we have 
Let us now prove that C'^^{f) = 0, for 4j + 1 < d. 

To do that it is convenient to introduce the following integral, for u > 0, w > 0, 
JK^fiu, v)=i J''-'^!^ fiz)dz. (4.21) 

We have easily 

Jk,ufiu,v)^^-^-^Jo,ufiu,v). (4.22) 
And using the residue theorem, we get 

A.fiu, v) = [i-ir+'fiv + zV^) + f{v - i^/u)) . (4.23) 
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From (|4.22p and (|4.23p we can compute Jk,ijf{u,v). 

To prove that c'^^{f) = for d > 4j + 1, we shall prove that each term in the 
sum (I4.20p vanishes, after integration in z and ^. 
Suppose first that j + 1 < fc < 2j + 1. We have 

Ql\x,P{x) - z,0 = 5]i?.,^(x)(P(x) - zYC- 

Hence 

2(PM-.)qy(^PW-..o 

IT L{z;x,^)''+^ 
is a sum of integrals like 



By integration by parts in z we have 



So, we can assume that 1^ = 0. But we have 

(-ly 

So we have I^{g){x,C) = near f = 0. Now we remark that for 

^ < 2j + 1 and d > 4j + 1 we have I < f + 1, hence ^ i-> I^{g){x,S,) is 
integrable and, using the analytic dilation argument already used for j — 0, we 
get 4{g)ix,^) = 0, hence 

2iP{x)-z)Ql\x,P{x)~z,0 



Now, assume that 2j + 2<k<3j. Using Lemma l3T2l we have 

Ql\x,P{x)^z,0^ J2 R.^^ix){P{x)-zrC- 

i^+|7|>2(2fe-l-2j) 

As above, we integrate by parts in z to have the possibility to put v at and 
then we use to decrease the order of the singularity in ^ as far as possible 

(integrability near ^ = 0) of ^l+i dz. We conclude by the analytic dilation 
argument. □ 

So, we have proven that in odd dimension d, c'^^ {f) = if 2j < 
We conjecture that the next following terms are not 0; more precisely we claim: 
Conjecture: For every j G N, j > 1, there exists / satisfying (|4.14p such that 
we have we have 

C^r" (/) ^ 0' and C^f (/) ^ (4.25) 



9 



In the following sections we shall check this conjecture for d = 1,3 and we 
shall compute analytic formula for cj^^\f) and c2\f)- Unfortunately, these 
analytic expressions have many terms and it is not obvious that (/) 7^ for 
d — 5,7, for every elliptic polynomial P. We shall see that this is true for convex 
polynomials for d — 7 and satisfying a technical condition if d = 5. Moreover 
we get, using numerical computations for particular non-convex polynomials P, 
that ^0. 

As we shall see in the next section, the property Cj^^ (/) 7^ gives easily a lower 
bounds on the density of eigenvalues. 

Remark 4.3 Following Remark \3.4\ we can extend our results to polyhomoge- 
neous polynomials P = P„i + Pm-i + ■ ■ ■ Pi + Pq- To follow the dependence in 
the coefficients, we note C'2j{,f,P) the coefficient C'2j{f) with polynomial P. 
In particular we have C'^^ {J , P^^^) — 0. /or d > 4j + 1 and for every e small 
enough. Assume now that d ~ 4jo — 3 or d — Ajo — 1, jo > 1. Then using a 
Taylor expansion in e, computed for e — h^/'--"^^^\ we get 

C,,{f,P^^'""'^''^)^Y.^un''^"'+'\ (4.26) 

*;>0 

in particular 70 = c'^^^(f,Pm) which is supposed to be not 0, as we have ex- 
plained before. 

5 Estimate the density of eigenvalues 

First of all let us remark that the nonlinear spectrum Sp[L] of L is included in 
the two quarters {z e C, 5R(z) > 0, ±3?(z) > 0}. 

On one side, it is easy to see that if A e K. and L{X)u = then u = Q. On the 
other side, if 3fJ(A) < and L(X)u = 0, computing Q{{L{X)u,u)) we conclude 
that u = 0. 

Let us denote by Nn{R) = #{z G Sp[L]; \z\ < R} and N{R) Nn=iiR). 

Proposition 5.1 For every real fi, jj, > d{m + l)/m, there exists C'k > such 
that 

NniR) < CkR^h-'^, Vi? > 1, Vft e]0, 1]. (5.27) 

V C'2j\f) with d > 2j, then for every r > 0, e > there exists c^,r > 
such that 

Nn{rh~') > c,^rh'\ Vfte]0,l], (5.28) 

where S = d — 2j . Moreover if j = (d even) then the estimates is valid with 
e = 0. So that, in even dimension, for every R> 0, Nti{R) behaves like h~'^. 

Proof. The proof of (|5.27p is a direct consequence of Weyl-Ky-Fan inequality 

m- 

We first remark that for every e > there exists R^ > such that if 
-7r/2-e < argz < 7r/2 + e, \^J.\ > Re 
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then we have 

\t + z\^>[l-e){t' + \z\^). 

Let us choose /(A) = (A + t)~^ with k large enough [^i > d{m + l)/m) and 
t > 0. We apply (|4.16p to get the following inequalities 

zeSp[L] zeSp[L] zeSp[L] 

But for every /i, ^i, large enough, such that /i — /ii is large enough, we have 

We choose now R — rh~^ to get 

Nnirh-') > ^ (1 + Ifiir" > c.^rh-' 

□ 

The above results concern the semi-classical regime. Now we give estimates 
for fi = 1 and high energy regime 

Corollary 5.2 For R /* +oo we have 

N{R) = C'(i?'i('"+i)/™). 

U C^if^if) 7^ with d—2j > 0, then for every £ > there exits > such that 
If j = 0, i/ie estimate is true with £ = and Cq > 0. 

6 1-d and 3-g? cases 

In this section we prove the following result. 

Theorem 6.1 For d = 1,3, there exits f satisfying Ii4.14\ l such that for every 
TO > 2, we have c'^\f) ^ 0. More precisely, we have 

&^\f) = -l-jj(^){P[x))P'{xfdx (6.29) 

Cf)(/) = I f'iPixmPix)\'dx (6.30) 

487r J^3 

We can choose /(A) = (A + i)^^ with fi > d{m + 1)/to and t > 0. 
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Proof. We compute with the exphcite form we got before for K2- We have, for 
d= 1,3, 

where H2 — 1 K2d^. But we have 

H2 = {P- zf-^Ap[b3 - 264a) +(P- zf-'^lWPf (63 - 264 - 2644) , 
hence 

^ (/) = ^ (^64,1 - ^63.0 + ^^4) ^ /(^^ {P{x))P'{xfdx. (6.31) 



c: 



if) - / f''Hp{x))p'{xrdx 
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and 



1 



487r ./„3 



(6.32) 



Ci'\f) = (2^(464,1 - 254 - 263) J^^ fiP{xmP\'{x)dx (6.33) 



(6.34) 



We have seen that for d odd, c? > 5, c'^\f) ~ 0. So we have to compute 
Cf^(/) for d = 5,7. 



7 5-(i and 7-d cases 



We have to compute in more details the term K4 from (j3.1ip . RecaU that we 
have 



Cf\f) = -2l 6{P-z){ I Ki{z-x,i)di\f{z)dzdx 



(7.35) 



We have to compute the foUowing three integrals, depending on a; G M*^ and 

ze C. 



(1) 



r(0;/3) 



1 



7 9; 



1 



i^:^ = r(a,o) 



. + {P{x) - zf « V lep + (^^(2:) - zf 
. I^P + (F(x) - z)' 



di;\ 



,dl{\iY)d<-,K2di- \a\ = 2, (7.36) 



r(3) 



r(0;/3) 



1 



+ {P(x) - z) 



(7.37) 
(7.38) 
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Using the new variable r] such that ^ = {P — z)ri (plus an analytic extension), 
we get 

a(/3) 7^ only when P — • • • , Pd) is such that /3j = 4 and /3fc = for A; 7^ j 

or /3j = = 2, j 7^ A; and /3£ = if € 7^ j, £7^/0. 

In the first case a(/3) = ai and in the second case a(/3) = 02 where 



96 (1 + \ri\^r 



It is convenient to introduce the following notations. 

_ f dr] _ f nfdT) _ f rifrifdri 

' ~ L (1 + w)^ ' " L (1 + wy ' ^'''^ " L (1 + • ^'^^ 

where j, fc, £ G N arc such that the integrals arc finite. Of course these integrals 
can be computed with the Euler beta and gamma special functions (see appendix 
for more explicite expressions). 

So we have a\ = 1^6,2 — l^s,! + ^'2 = ^'e.i.i- 

Using integration by parts, in x or in ^, we get the following formulas 

C4(/)= / Ciif;x)dx (7.43) 

where 



and 



C4f; x) = C4,i(/; x) + CiMf; x) + C^Mf; (7-44) 
C4,i{f;x) = 2(f{P-z)Y.[dS{P-zfl^^^)f{z)dz, (7.45) 

^ |/3|=4 

2 j2 ^("''^^i ^1^''^^" {^(i) ^"^^) 



\a\=2 



C^,^{f;x) = 2£{P-z)dSL{z)lff{z)dz. (7.47) 

Now we compute each term. After elementary but tedious computations we get 
the following results, using the notations: 



ai = ^66,2 - ^&5,i + -^^'4 (7.49) 
a2 = be,i,i (7.50) 
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For d=5 , we have 



where 
El 
E2 
E3 



G2 



G3 



l<j<5 j<k 

l<j<5 j<k 



(7.51) 



C4,2(/;x) 2--' 4 

cM^) = ^f{p)-^np)+^f^'\p)- 



12' 
12' 



(APf (64-265-265,1 + 466,1), 
|VP|2Ap(64 + 1266 - IO65 - 265,1 + 1666,1 - 1667,2), 
|VP|''(2066 - 665 - I667 + 1667,1 - 1266,i), 
(AP)2(l266,i - 265) - 1667,2^(^1^)' 



(7.52) 
(7.53) 

(7.54) 
(7.55) 
(7.56) 
(7.57) 



-3267,1.1 ^(a|p)(a^p) - 1667,1,1 5^ (5^,p)^ 

j=^k j=ik 

API VPp (2466,1 - 1667,1 - 465 + 465) - 3267,2 ^(ajP2)(ajP)' 

3 

-3267,1,1 J2i9jP)idkPf - 6467,1,1 J2i9lk)idjP){dkP), (7.58) 

j^k j^k 

|VP|^(l266,i - 1667,1 - 265 + 466) - 1667,2 ^(a,P)^ 

j 

-4867,1,1 ^(a,p)2(afcP)2. (7.59) 

Finally, we get 

Ciif;x) = Ao{x)f{Pix))+A,{x)f'iPix)) + A,ix)nPix))+A3ix)f^^\P{x)), 

(7.60) 
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Ao{x) = -20ai^dfP-20a2j29jdlP, (7-61) 

A,{x) = 8(ai-6r,2)I](5|P)2 + 4(a2-26r,i)^(5|,,P)2 (7.62) 
{APf 



(&4 - 465 - 265,1 + 1666,i) - 1667,1,1 Y,i9jPMP), 

A2ix) = J ^ (^1465-1666-64 + 265,1-4066,1 + 3267,1] (7.63) 

+867,2 Y.{d]PMPf + 867,1,1 J2idjP)idkPr 

+1667,1,1 ^(a^,p)(a,p)(9,p), 

^3(0;) = ^^(2466-865-1667) -^67,2 X](a,P)" 

j 

-467,1,1 ^(a,p)2(afeP)2. (7.64) 

Using explicit computations (see Appendix), we get 

j i<fc 

_i[l(AP)2 - ^ y (a?p)(a2p), 

160^ ' 240^^ ^ A fe 

A2(x) = ^VPrAP+^^(a|P)(a,P)2 (7.67) 

j 

+^ T.idlp)idkPf + ^ Y.^di,p){d,p){d,p), 
Mix) = -^ivpr-^E(^.-^)^-^E(^.-^)^(^^^)M7.68) 



For d =7 1 we have 

C4f; x) = Ao{x)f{P{x)) + ^i/'(P(a;)), (7.69) 
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IVPrAP/ \ 
Mx) = ^ (14^5 - 1666 -bi + 265,1 - 4066,1 + 3267,1) (7.70) 

+16br,2j2i9]P)id,Pf + 1667,1,1 ^(a|P)(afeP)2 

+3267,1,1 ^(a^feP)(5,p)(afeP), 

Mx) = (2466 - 865 -I667) -867,2 ^(SjP)^ 

3 

-2467,1,1 ^(a,P)2(9feP)2. (7.71) 
Using explicit computations as for = 5 (see Appendix), we get 

^O(-) = ^^T.^9-3P)id3Pf + ^^T.^d^P){duP? 

+ ^T.(9lkPMPmP), (7.72) 

Mx) = -^ivpr-^E(^.-^)'-|^E(^.^)'(^^^)M7.73) 

3 3i^k 

We should like to use these formulas with /(A) = (A + i)"**, /U > d(jn + l)/m 
and t > 0, to prove that cf''{j) = /g^ C4(/; x)dx ^ (rf = 5, 7). 

With / like above, we can see easily that C^f^ (/) ^ for the following polyno- 
mials: 

|T| P{x) = E a^ajf , Q!j > 0, 1 < j < d, d = 5, 7. 

i<j<d 

[2] P(a:;) = aj^kXjXk, is a positive-definite quadratic form, d = 5, 7 

l<j,fe<d 

[3] d = 7 and P is convex. 

[4] d = 5, P is convex and satisfies the inequalities 

E dplP < 2j2 9tP (7.74) 

l<j<fc<5 l<j<5 

j:{di,py < im y: {d^pf (7.75) 

j¥fc l<j<5 

For non-convex polynomials, wc can check that C^'^\f) =/= 0, d — 5,7, for many 
examples with numerical computations, supporting our conjecture that for every 
elliptic polynomial P, C^^\j) ^ 0, if d = 5, 7(see Appendix). 

For d = 9, 11, it seems difficult to compute CQ^\f) by hand. We need more 
help from symbolic and numerical computations to check our conjecture. 
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8 appendix 

8.1 Formulas for bj^k,e 

We assume d> 3 and 2j — q > 1. We have 



where 



/■+°° r« 1 g- 

Jo (TT7^'^=2^^^ 

B{x,y)= [\^-\l-t)y-' 
Jo 



-,J 



r{x + y) 



So computing in polar coordinates, 



bj{d) = 



TT'^/^rp- - d/2) 

m 



Now, by elementary computations we get easily 



(8.76) 



bjAd) = -{bj.i{d)-bj{d)) 



i?(5/2,j-^)6,(d-l) 
lB{3,j-^Md-2). 



(8.77) 

(8.78) 
(8.79) 



8.2 Numerical computations for C4(/) 

The following computations has been performed by Guy Moebs, Research En- 
gineer, Laboratoire Jean-Leray, CNRS-University of Nantes. 

The method used to compute multi-dimensional integrals is Montc-Carlo, 
with a cut-off of the domain to reduce it in a bounded domain fitting with the 
behaviour of the polynomial P. 

In each example, 100 simulations are computed with at least 10® events. 

d 

Example. 1 d= 5, polynome P(x) = ^^xj^ + axi^X2^ 



a 


C4{.f) 


7 


1 428 


10 


1 515 


100 


9 237 


1000 


235 115 
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d 

Example.2 d=l, polynome P(x) = ^xj^ + q;xi^X2^ + /3x3^X4^ 

j=i 



a 


13 


C4(/) 


7 


7 


409 


7 


10 


423 


7 


100 


1 806 


7 


1000 


39 646 


10 


10 


434 


10 


100 


1 705 


10 


1000 


36 724 


100 


100 


1 755 


100 


1000 


19 587 


1000 


1000 


18 270 



d 

Example. 3 d = 5, polynome P(x) = ^^xj^ + axi^X2^ + /3x3^X4^ 

j=i 



{a, (3) 


C4(/) 


(100, 10) 


11 732 
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